Time-resolved photoemission of correlated electrons driven out of equilibrium 
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We describe the temporal evolution of the time-resolved photoemission response of the spinless 
Falicov-Kimball model driven out of equilibrium by strong applied fields. The model is one of the 
few possessing a metal-insulator transition and admitting an exact solution in the time domain. The 
nonequilibrium dynamics, evaluated using an extension of dynamical mean-field theory, show how 
the driven system differs from two common viewpoints - a quasiequilibrium system at an elevated 
effective temperature (the "hot" electron model) or a rapid interaction quench ("melting" of the 
Mott gap) - due to the rearrangement of electronic states and redistribution of spectral weight. The 
results demonstrate the inherent trade-off between energy and time resolution accompanying the 
finite width probe pulses, characteristic of those employed in pump-probe time-domain experiments, 
which can be used to focus attention on different aspects of the dynamics near the transition. 

PACS numbers: 71.10.Fd, 78.47.J-, 79.60.-i, 03.75.-b 
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I. INTRODUCTION 

Recently, pump-probe techniques, successfully em- 
ployed in optical reflectivity studies ?i have been used 
to extend photoemission spectroscopy 2 - (PES) to the 
time domain (time-resolved PES or tr-PES) in the fem- 
tosecond and attosecond regimes.— ~— These advances 
open the possibility of observing and controlling dynam- 
ics on time-scales relevant to correlated electronic pro- 
cesses^— specifically in optical pump-probe radio fre- 
quency (rf) cold atom and table-top laser experiments, 
as well as using free-electron laser facilities, such as the 
Free-electron LASer in Hamburg (FLASH) and the Linac 
Coherent Light Source (LCLS), to conduct pump- probe 
extreme ultra-violet or x-ray photoemission and scatter- 
ing studies in the time domain. 

Most pump-probe experiments have been interpreted 
in terms of "hot" electrons, effectively equilibrated at 
highly elevated temperatures^—^ While this approach 
can capture a redistribution of spectral intensity through 
the change in the Fermi distribution function and ther- 
mal rearrangement of electronic states that experiments 
observe on picosecond time scales, it can not account 
properly for the out-of-equilibrium rearrangement of ac- 
cessible electronic states nor capture the nonequilibrium 
redistribution of spectral weight that accompanies pump 
pulses with the high excitation densities needed to drive 
phase transitions or excite collective modes characteristic 
of correlated electron systems on ultrashort time scales 
in the femtosecond or attosecond regimes. 

As a test case, we have chosen to study the effect 
of strong driving fields on a simple model system for 
which the hot electron model definitively breaks down 
and where the effect of the driving field does not mimic 
an interaction quench 1 * 3 - - — or "melting" of the Mott gap. 
We examine the temporal evolution of the tr-PES re- 



sponse for the spinless Falicov-Kimball model at half- 
filling, driven by a large amplitude, dc electric field to- 
ward a nonequilibrium steady-state , d 8 The Falicov- 
Kimball model is one of the simplest correlated electron 
models and it has a Mott-Hubbard metal-insulator tran- 
sition (MIT) at half-filling; to this point, it is the only 
model for which an exact nonequilibrium impurity solu- 
tion has been developed in time-dependent fields with a 
time range long enough to evaluate tr-PES. In particu- 
lar, the temperature invariance of the equilibrium den- 
sity of states (DOS) for this model makes a comparison 
to hot electrons at long time delays relatively straightfor- 
ward. We find that the spectral intensity develops regular 
Bloch-like oscillations for weak metallic correlations that 
become sharply damped approaching and crossing the 
MIT. The results elucidate the out-of-equilibrium behav- 
ior of a simple correlated electron system observed using 
tr-PES as a probe. 



II. METHOD 

We determine the real time dynamics for the model 
on the hypercubic lattice in infinite dimensions (d — 
oo) using nonequilibrium dynamical mean-field theory 
(DMFT^iiLiS This method yields the double-time 
contour-order Green's function (GF) G (t,t') within the 
Kadanoff-Baym-Keldysh formalismj 19 i 20 The system be- 
gins in thermal equilibrium at time t = i m j n and temper- 
ature T before an electric field, applied at t = 0, breaks 
time-translation invariance. The system evolves under 
the influence of this field to a maximal time t max . This 
defines the Keldysh contour C used in the formalism. The 
contour-ordered GF encodes both the retarded GF, de- 
termining the equilibrium as well as nonequilibrium ar- 
rangement of states, and the lesser GF, specifying the dis- 
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tribution of electrons among these states, together with 
other physically relevant GFs. In particular, the lesser 
GF, related to the PES response, is given by the Keldysh 
contour-ordered quantity, 



GfAt,t') = iTr 



exp 



-PH K 



c](t') Ci (t)/Z eq 



in the Heisenberg picture, where t lies on the upper real 
time branch and t' on the lower real time branch of the 
Keldysh contour and "eq" denotes equilibrium (i < 0) 
quantities at the initial temperature T. 
The equilibrium Hamiltonian is given by 
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describing the hopping, t* , of conduction electrons be- 
tween lattice sites with a filling controlled by the chemi- 
cal potential, fj,, that experience an on-site interaction, U, 
with another species of localized electrons with an occu- 
pation Wi. U c = v2 1* is the critical interaction strength 
for the MIT at half-filling. Throughout this work, the 
energy unit is taken to be t* . As an aid to the reader 
in understanding the relevant time and field scales for 
this paper, consider those set by the hopping integral 
t* = 250 meV and hypercubic lattice spacing a — 3 A; 
the corresponding unit of time is ~ 16 fs and that for 
the dc driving field E is ~ 13 mV/A. The nonequilib- 
rium DMFT formalism proceeds in essentially the same 
manner as the iterative approach applied in equilibrium 21 
where all quantities now have two time indices.— 

The driving term is modeled by a spatially uniform 
constant dc electric field along the (1, 1,1,.. .) hypercubic 
body diagonal high symmetry direction to simplify evalu- 
ation of the noninteracting GFJ£ The spatially uniform 
vector potential, in a gauge with zero scalar potential 
(Hamiltonian gauge), associated with this driving field 
varies linearly in time and enters through a Peierls' sub- 
stitution.^ 

We determine the real frequency spectral intensity us- 
ing a finite width probe pulse that samples the real 
time dynamics of the driven nonequilibrium system. The 
probe pulse envelope in a tr-PES experiment can be well 
approximated by a Gaussian waveform 



s(t) 
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where t Q measures the time delay with respect to the ap- 
plication of the driving field and er measures its effective 
temporal width. The tr-PES response function is then a 
probe pulse weighted relative time Fourier transform of 
the lesser GFM^^. 



III. RESULTS AND DISCUSSION 

Results are shown in Fig. Q] for metallic systems with 
weak correlations U < U c [(a) - (c) U = 0.125 and (d) 



.. (a) (c) 

-1 ^^^^^^^^^^^^^B ^^^^^^^^^^^^^^ 

an'i -, ' , ! 1 m ' 1 i 

a 2 A A A . 



w 1 



-0.4 
J 0.3 

1 0.2 

: : 0.1 
-0.0 



"0 10 20 30 40 50 10 20 30 40 50 

Time Delay (1/t*) 

FIG. 1. (Color online.) Tr-PES intensity (falsecolor or 
grayscale) as a function of photoelectron energy and time de- 
lay for weakly correlated metallic systems U < U c [(a) - (c) 
U = 0.125 and (d) U = 0.25]. A probe pulse of characteristic 
width a — 1.0 samples the nonequilibrium dynamics of the 
system driven by fields of strength (a) E = 0.5, (b) E = 1.0, 
and [(c) and (d)] E — 2.0. Bloch oscillations with a period 
proportional to 1/E develop almost immediately following ap- 
plication of the field at t = and an additional amplitude 
modulation with a period proportional to 1/(7 appears for 
the strongest driving fields [(c) and (d)]. 



U = 0.25] driven by applied fields with different strengths 
E. Each panel shows the spectral intensity (falsecolor 
or grayscale), as a function of binding energy and time 
delay (obtained from data generated with a discretiza- 
tion in time equal to 0.1). For the lowest field strengths 
[Fig. (Ha), E = 0.5, and Fig. \V[b), E = 1.0] the spec- 
tral intensity develops regular Bloch oscillations, with 
a period proportional to 1/E, damped by correlations. 
For sufficiently large fields [Fig. Hie), U — 0.125, and 
Fig. Did), U = 0.25, with E = 2.0], an additional am- 
plitude mode, characterized by "beats" in the spectral 
intensity, appears with a period proportional to 1/U. 

Figure [5] shows results for interaction strengths ap- 
proaching and crossing the MIT. Oscillations associated 
with the driving field are still apparent for metallic sys- 
tems U <U C [Fig.rjfa), U = 0.5, and Fig.^b), U = 1.0], 
but the amplitude mode and increased damping lead to 
a rather irregular temporal evolution. For U = 1.5, just 
above U c , the damping is severe enough to suppress os- 
cillations for all but the shortest time delays. Further 
increase in the interaction strength exacerbates these ef- 
fects. The behavior of these oscillations is the tr-PES 
analog of that found for the instantaneous current re- 
sponse, evaluated from the equal-time lesser GFJ^ 

The results shown in Figs, [^c) and [UJd) naively sug- 
gest, at least for U > U c , that the dc field drives the 
system toward a nonequilibrium steady-state in which 
the Mott gap has melted or the interaction strength has 
been quenched to a smaller value U < U c , resulting in 
significant spectral weight at and above the equilibrium 
Fermi level (Energy = 0) for long time delays. However, 
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FIG. 2. (Color online.) Tr-PES intensity (falsecolor or 
grayscale) for various correlations in both metallic U < U c 
[(a) U = 0.5 and (b) U = 1.0] and insulating U > U c [(c) 
U = 1.5 and (d) U = 2.0] regimes driven by fields of strength 
E — 2.0 and sampled using a probe pulse of width a — 1.0. 
Bloch oscillations become increasingly damped approaching 
and crossing the MIT. 



these observations are merely artifacts of the trade-off 
between energy and time resolution associated with the 
relatively narrow Gaussian probe pulses used in Fig. [5] 
to highlight the temporal dynamics. In the transient re- 
sponse regime there is no time translation invariance and 
the probe width affects both the temporal and energy res- 
olution. However, the conventional Fourier uncertainty 
relations would be recovered in the steady-state regime 
at long timesi 14 ! 24 

Figure [3] shows results similar to those of Fig. (2J but 
for wider probe pulses. Note that the increase in width 
leads to a suppression of temporal dynamics, except at 
the shortest time delays, and sharper spectral features 
as a function of energy. For the weakly correlated metal 
[Figs. [3Ja) and EJe), U — 0.5], the increased width sup- 
presses the regular Bloch oscillations. The spectral in- 
tensity approaches the steady-state at long times includ- 
ing the redistribution of weight into faint high- and low- 
energy satellites at ~ ±E associated with the rungs of the 
Wannier-Stark (WS) ladder^ The WS ladder describes 
the rearrangement of electronic states into periodic reso- 
nances in energy within a conventional metal or weakly 
interacting system due to the application of a strong driv- 
ing field. The spectral intensity for the strongly corre- 
lated metal [Figs. |3Jb) and EJf), U — 1.0] behaves simi- 
larly with the appearance of additional WS satellites and 
a suppression of weight near zero energy attributable to 
the rearrangement of accessible electronic states 1 ^ and 
not an ad hoc change to the interaction strength. 

For U ~ U c [Figs. |5] (c) and [31(g)] there is still a signifi- 
cant redistribution of weight across the equilibrium Mott 
gap at short time delays before relaxing and partially re- 
covering at longer times. Further increase in the probe 
pulse width (a = 4.0) suppresses the temporal evolution 
but does reveal a modified real frequency structure. Fig- 
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FIG. 3. (Color online.) Tr-PES intensities (falsecolor or 
grayscale) from the metallic to insulating regimes for different 
probe pulse widths. Together with increased energy resolu- 
tion, wider probe pulses lead to qualitative changes in the 
temporal evolution of the PES response and to a sharpening 
of spectral features as a function of energy. The parameters 
are [(a) and (e)] U = 0.5, [(b) and (f)] U = 1.0, [(c) and 
(g)] U = 1.5, and [(d) and (h)] U = 2.0 for a = 2.0 and 4.0, 
respectively, all driven by fields of strength E = 2.0. 



uresO (d) and (h) show results for U = 2.0. In this 
case, spectral weight is distributed between two main fea- 
tures and the temporal evolution appears only through 
a modulation of the spectral intensity within these fea- 
tures. The ability to track changes in the spectral inten- 
sity persists only for the shortest time delays after ap- 
plying the driving field and the spectra quickly approach 
that characteristic of the nonequilibrium steady-state. In 
each of the cases presented in Fig. [21 the tr-PES response 
clearly is not indicative of an interaction quench with 
subsequent melting of the equilibrium Mott gap. 

Finally, we compare the response in equilibrium (taken 
at time delay t a = —1) to the response approaching the 
regime of the nonequilibrium steady-state (taken at time 
delay t Q = 15). Figure 2] shows this comparison for both 
a metallic [Fig. Q|a), U = 0.5] and insulating [Fig.Hfb), 
U = 1.5] system probed by pulses of widths a = 2 and 
4. The response in equilibrium (black curves) essentially 
matches the equilibrium DOS multiplied by the Fermi 
distribution function and convolved with an energy reso- 
lution function accounting for the finite temporal width 
of the probe pulse. The energy resolution improves with 
wider probe pulses, manifest in a sharpening of spectral 
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FIG. 4. (Color online.) Equilibrium and nonequilibrium 
steady-state response for (a) U — 0.5 and (b) U = 1.5 with 
a driving field strength E = 2. These results correspond to 
time delay cuts at t a = —1 and t = 15 from the results 
shown in Figs. EJa), [3jc), [3je), and Etg), sampled with two 
different probe pulses of width a — 2.0 and a — 4.0, respec- 
tively. The nonequilibrium steady-state response (time delay 
t = 15) shows high and low energy Wannier-Stark satellites 
(highlighted by arrows) and overall line shape, including a 
sharpening of spectral features, incompatible with both the 
simple hot electron view and melting of the Mott gap or an 
interaction quench with applied field. 

features between a = 2 and 4. 

The response in the nonequilibrium steady-state (red 
or gray curves) shows qualitative differences to those 
in equilibrium. Note that the equilibrium DOS in the 
Falicov-Kimball model at half-filling is symmetric with 
respect to the Fermi level and temperature independent; 
therefore, at a higher effective temperature, the spectral 
weight should be redistributed to at most one higher- 
energy feature above the equilibrium Fermi level. How- 
ever, in these cases both high and low energy WS satel- 
lites are found in the response, indicated by arrows, and 
the features are far sharper than those in equilibrium. 
For the insulator shown in Fig. HJb), there is even a 
suppression of spectral weight at the Fermi level. The 
satellites are more pronounced for weak correlations and 
wider probe pulses [see Fig. [3] and compare Figs. Ufa) 



and HJb)]; however, they are present in the response for 
all cases, highlighted here by systems on both sides of the 
MIT. Overall, this behavior precludes a simple quasiequi- 
librium description of the out-of-equilibrium response of 
the system in terms of an elevated effective temperature 
- the "hot" electron model. 



IV. CONCLUSIONS 

The current model captures the formation of damped 
Bloch oscillations for weakly correlated metals. The os- 
cillations, typically suppressed in real materials due to 
scattering from phonons and impurities, not included in 
this model, are simply damped here by electron corre- 
lations. Experiments corresponding to the conditions 
represented in these simple model calculations poten- 
tially could be performed in ultracold atomic systems 
by generalizing equilibrium rf techniques^ to nonequi- 
librium situations. The experiment would involve mix- 
tures of light fermions with heavy fermions or heavy 
bosons at low temperature (Li 6 -K 40 mixtures for the for- 
mer or Li 6 -Cs 133 or K 40 -Rb 87 mixtures for the latter), 
applying a driving field to generate the Bloch oscilla- 
tions. This "field" could be gravity or a detuning of 
the counter-propagating lasers which "pulls" the optical 
lattice through the atomic clouds. The rf pulses would 
need to have a duration which is short enough in time 
to observe the Bloch oscillations in the time domain af- 
ter a time-of-flight image. This type of experiment may 
be cleaner than those performed on conventional con- 
densed matter systems because the driving fields will not 
interfere with the time-of-flight detection used to probe 
these systems. For conventional condensed matter sys- 
tems, using ultrafast probe pulses to determine the tr- 
PES response of a system driven by a strong electric field 
toward a nonequilibrium steady-state may be challeng- 
ing to replicate experimentally. The ability to observe 
short-time behavior (on the scale of femtoseconds) opens 
the possibility of observing oscillations before they be- 
come damped (something that may occur at picosecond 
time scales, especially in weakly correlated materials). It 
is conceivable that the WS ladder or Bloch oscillations 
could be seen within the duration of a wider pump exci- 
tation using a particularly sharp probe pulse. This could 
be true for FEL sources with exceptionally large ampli- 
tude pump pulses. However, this method would require a 
modification of existing synchronization techniques and 
improvements in probe temporal resolution to reach the 
necessary time scales. 
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